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ABSTRACT 


A  parametric  analysis  is  given  of  imaging  in  backscattered  light  in 
clear  ocean  water  characterized  by  temperature  and  salinity  inhomogeneities 
whose  scale  length  is  large  c  ..  pared  to  the  diameter  of  the  light  beam.  An 
expression  for  the  beam  pattern  of  a  finite  underwater  laser  beam  is  derived 
where  the  effects  of  source  coherence  on  the  imaging  properties  of  back- 
scattered  light  are  discussed.  On  the  basis  of  the  model  assumed  here,  it  is 
concluded  that  the  resolution  cannot  be  significantly  improved  by  redaction  of 
the  camera  threshold.  This  limitation  results  from  the  rapid  dem-.ase  of  the 
modulation  transfer  function  for  increasing  spatial  frequencies,  ’  ,nich  indi¬ 
cates  that  image  modulation  below  two  percent  falls  another  10  dl '  jt  only 

negligibly  higher  spatial  frequencies.  Thus,  a  10-dB  improver . t  m  the 

limiting  resolution  of  the  camera  will  not  result  in  significant!^  better  image 
quality  through  ocean  water. 
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I.  INTRODUCTION 


The  imaging  pioperties  of  collimated  laser  beams  propagating  through 
ocean  water  have  recently  attracted  considerable  attention.  Direct  or  indi¬ 
rect  viewing  is  necessa-v  to  carry  out  such  underwater  activities  as  target 
acquisition,  guidance  of  submarines,  and  the  seaich  for  metal  nodules  or  ma¬ 
rine  specimens.  Clear  ;cean  water  (i.  e.  ,  free  of  particulate  matter)  is  in¬ 
homogeneous  because  of  the  existence  of  random  variatiens  in  salinity  and 
temperature.  Both  quantities  affect  the  index  of  refraction  of  wa.'er;  conse¬ 
quently,  a  laser  beam  employed  to  illuminate  a  subsurface  object,  even  in 
clear  ocean  water,  may  be  subject  to  severe  degradation.  Whether  a  laser 
beam  is  to  be  used  for  ranging  or  illumination  (in  direct  or  in  TV  viewing)  or 
for  communication,  knowledge  of  its  behavior  is  essential.  The  complex 
interrelationships  among  natural  limitations  and  the  necessity  of  obtaining  a 
clear  image  impose  severe  restrictions  on  designers  of  underwater  equipment 
that  employ  coherent  light  beams.  In  this  paper,  mathematical  expressions 
are  presented  that  permit  estimation  of  the  performance  of  proposed  equip¬ 
ment  in  clear  ocean  water  under  various  operational  conditions. 

In  particular,  this  paper  presents  a  parametric  analysis  of  underwater 
imaging  in  backscattered  light  in  ocean  waters  that  are  relatively  free  of  par¬ 
ticulate  matter.  In  this  case,  light  is  scattered  by  refractive-index  variations 
in  the  water  resulting  from  large-scale  (some  tens  of  centimeters)  thermal 
and  saline  variations  (about  their  respective  means).  Ocean  waters  of  this 

type,  which  'nave  long  attenuation  length-*  (~8  in)  in  the  blue-green  region  of 

1  2 

the  spectrum,  occur  in  several  geographical  locations:  ’  e.g.,  the  Caribbean 

Sea,  the  Pacific  North  Equatorial  Current,  the  Pacific  Countercurrent,  the 
Pacific  Equatorial  Divergence,  the  Pacific  South  Equatorial  Current,  and  the 
Sargasso  Sea. 

The  mutual-coherence  function,  defined  as  the  cross -corre1ation  func¬ 
tion  of  the  complex  field  in  a  direction  transverse  to  the  direction  of  propaga¬ 
tion,  describes  the  loss  of  coherence  of  an  initially  coherent  wave  propagating 
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in  an  inhomogeneous  medium.  It  follows  that  the  mutual-coherence  function 
is  impo--_.it  lor  a  number  of  practical  applications.  It  determines  the  limit¬ 
ing  resolution  obtainable  in  forming  an  image  through  an  inhomogeneous  me¬ 
dium,  the  x’nean  intensity  distribution  from  an  initially  coherent  wave  emanat¬ 
ing  from  a  finite  transmitting  aperture,  the  limiting  holographic  resolution  in 
an  inhomogeneous  medium,  and  the  signal-to-noise  ratio  in  a  heterodyne  detec¬ 
tion  system. 

Section  II  presents  a  parametric  analysis  of  visibility  in  clear  ocean 
water.  The  effects  jf  source  coherence  on  the  imaging  properties  of  back- 
scattered  light  are  discussed.  In  addition,  an  expression  for  the  beam  pattern 
of  an  underwater  laser  beam  is  presented,  and  the  quantitative  relationships 
of  image  resolution,  beam  .-pread,  and  temperature  and  salinity  fluctuations 
are  di«  cussed.  Finally,  an  expression  for  the  limiting  resolution  in  clear 
ocean  water  is  presented. 
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II.  COHERENCE  EFFECTS 


A.  PARAMETRIC  ANALYSIS  OF  VISIBILITY  IN  CLEAR 

OCEAN  WATER 

This  paper  treats  time -averaged  quantities  in  which,  for  a  homogeneous 
isotropic  turbulent  medium,  the  mutual  coherence  function  is  real;  this  func¬ 
tion  is  referred  to  here  as  the  modulation  transfer  function  (MTF).  Further¬ 
more,  we  invoke  the  ergodic  hypothesis,  in  which  the  time  average  and  the 
ensemble  average  are  assumed  equal.  The  calculations  that  follow  pertain  to 
the  case  of  large-scale  (in  comparison  with  laser  beam  diameter)  thermal  and 
saline  variations  about  their  respective  means. 

Yura^  has  given  a  quantitative  analysis  of  small-angle  (£  1  mrad)  scat¬ 
tering  of  light  in  ocean  water  by  suspended  biological  particles  with  an  index 
of  refraction  close  to  that  of  water  and  large-scale  (in  comparison  with  the 
laser-beam  diameter,  £  1  cm)  index-of-refraction  variations  due  to  thermal 
and  saline  inhomogeneities.  In  particular,  he  has  calculated,  for  plane  waves, 
the  MTF  for  these  two  scattering  mechanisms  and  has  shown  that  (for  values 
of  the  transverse  distance  p  less  than  the  size  of  the  large-scale  index-of- 
refraction  fluctuations  but  larger  than  the  size  of  the  suspended  biological 
particles)  the  MTF's  due  to  the  two  mechanisms  have  different  functional 

dependences  on  transverse  distance.  In  particular,  for  the  case  of  large- 

3 

scale  index-of-refraction  variations,  he  finds  that 


M(p,  R)  =  exp 


k2p2<An2>R 

4a 


(1) 


where  ^An"^  is  the  rms  index-of-refraction  variation,  a  is  the  character istic 
scale  length  of  the  inhomogeneity  (assumed  greater  than  the  laser  beam  diam¬ 
eter),  p  is  the  transverse  distance  at  propagation  distance  R,  k  is  the  optical 
wave  number  (2rm/\),  and  angular  brackets  denote  the  ensemble  average. 
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For  calculation  of  such  quantities  as  resolution  and  beam  spreading, 
the  MTF  of  a  point  source  (i.e.,  spherical  wave)  must  be  known.  In  this 
case,  the  spherical-wave  MTF  is  obtained  from  the  plane -wave  MTF  by 
replacement  of  p^  by  (pu)^  du,  which  equals  p^/3.^ 

Ms(p,R)  ^  exp  (2) 

2 

The  mean-square  index  variation  (An  >  is  given  by 

<An2>=  (| |)2<AT2>  +  (|§)2<  AS2>  (3) 

2  2 

where  (AT  >  and  <AS  >  are  the  mean-square  temperature  and  salinity  vari¬ 
ations  (about  their  respective  means).  In  deriving  Eq.  (3),  we  have  assumed 

that,  to  a  first  approximation,  temperature  and  salinity  are  independent 

2  1/2 

random  variables .  Typically,  <AT  >  '  is  a  few  hundredths  to  a  few  tenths 

5  2 

of  a  degree,  absolute;  however,  no  data  are  readily  available  on  <AS  >. 

Also,  the  characteristic  scale  length  a  for  temperature  fluctuations  is  on  the 

5 

order  of  tens  of  centimeters. 

In  general,  the  incremental  index-of- refraction  variation  An  is  related 
to  the  incremental  density  variation  Ad  by 

■  (fX^H)^)  <4> 

3 

where  d  is  the  density  (  =  1.02  g/cm  for  ocean  water).  In  deriving  Eq.  (4), 
we  have  neglected  the  change  of  polarizability  of  the  constituent  molecules 
with  changes  in  density.  For  the  values  of  the  relative  temperature  and 
salinity  variations  considered  here,  order-oi-magnitude  estimates  of  the 
resulting  change  in  the  polarizability  show  that  this  neglect  is  justified. 
Assuming  that  tempeiature  and  saJimty  variations  are  independent  variables, 
we  may  write  Ad  -  (Ad)^,  +  (Ad)^,  whej  e 


V* 
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=  -a  AT 


(5) 


-4  o 

where  O'  is  the  coefficient  of  volume  expa  -.ion  (~2.  i  X  10  /  C  for  water  in 

the  range  0  to  33  °C),  and  (Ad)g  =  A(dg),  where  dg  is  the  saline  contribution  to 
the  density  (i.e.,  the  concentration  of  the  salt  times  its  density).  Hence,  on 
squaring  and  taking  an  ensemble  average,  we  obtain 

<a"2>  *  {  P<at2>]  +  to.)2]  <6> 

Since  no  data  are  readily  available  on  Ad,.,  we  consider  temperature  variations 

L  ‘ 

only  in  the  following.  In  this  case  we  find  that,  for  ocean  water  (n  -  1.34), 

<An2>  =  0.6  X  10’8<AT2>  (7) 

For  example,  for  10 <  <An2>  <10  ^ ,  we  find  from  Eq.  (7)  that  0.  13  °C 
2  1/2 

<<AT*>  <  0.41  °C.  (Since  saline  variations  have  been  neglected,  these 

rms  temperature  variations  are  to  be  regarded  as  an  upper  limit.) 

Figure  1  is  a  graph  of  the  spherical  wave  MTF  [i.e.,  Eq.  (2)]  for 
X  -  0.  5  pm  over  the  range  of  interest.  The  upper  abscissa  is  given  for  p  in 
centimeters,  while  the  lower  one  is  given  for  spatial  frequency  f  in  cycles/ 
milliradian  (f  =  p/2\).  Values  of  the  parameters  (An2^and  scale  length  a 
were  chosen  for  which  the  MTF  would  be  about  2  percent  at  about  1 . 0  to  10.0 
cycles/mrad  and  X  =  0.5  pm,  with  the  constraint  that  the  scale  length  a  have  a 
reasonable  value.  The  curves  are  shown  for  a  =  50  cm  and  (An  ^  -  10  and 


Receiving  systems  differ  greatly  in  their  resolution  thresholds.  Further¬ 
more,  these  thresholds  are  a  function  of  contrast^  and  level  of  illumination 

7 

reaching  the  receiver.  However,  even  a  few  tenths  of  a  watt  of  laser  illum¬ 
ination  should  be  adequate  to  realize  the  limit -ug  resolution  of  current  low- 
level  television  systems.  For  any  given  receiving  system  there  is  seme 
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10 1  10°  10' 
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Figure  1.  Spherical- Wave  MTF  as  Function  of  Transverse 
Separation  p  and  Spatial  Frequency  f  for  Various 
Va'ues  of  R  and  <An^> 
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number  mQ  (<1)  for  which,  if  the  MTF  is  less  than  m0,  resolution  is  not 
possible.  Defining  pmax  so  that,  for  a  fixed  propagation  distance  R, 

M(p  ^  Pmax)  ~  mo>  we  n°te  that  the  minimum  distance  resolvable  in  the 
object  plane  AXjj^  is  related  to  the  maximum  transverse  distance  in  the 
receiving  aperture  pmax  (assuming  pmax  <  receiver  diameter)  by 


Ax  .  = 


mm 


R 


kp 


max 


(8) 


The  contrast  degradation  caused  by  the  medium  is  described  by  Fig.  1,  in 

7  1 

which  the  2 -percent  modulation  level  is  marked  by  the  heavy  horizontal  line. 
From  Fig.  1  we  see  that,  for  the  case  of  large-scale  index  variations, 

M  5  2  percent  for  0.1  <  p  <  1.0  cm'. 

In  general,  to  determine  ihe  effective  resolution  through  such  a  medium, 
one  should  plot  the  limiting  resolution  as  a  function  of  apparent  contrast  for 
the  receiver  of  interest.  The  intersection  between  the  limiting  resolution  , 
curve  of  the  camera  and  the  MTF  of  the  medium  is  the  achievable  resolution. 
For  example,  at  7-percent  contrast  (about  3.5-percent  modulation)  and  high 
light  levels,  the  limiting  resolution  for  a  1-in.  image1  orthicon ,  S-10  photo- 
cathode,  f / 1.5  optics,  6-Mc  gain  bandwidth,  is  about  16  TV  lines/mm.  0,7 
With  a  1-m  focal  length,  this  corresponds  to  a  spatial  frequency  of  8  cycles/ 
mrad.  But  from  Fig.  1  we  see  that,  except  for  (An  >  =  10  at  a  range  of 
10  m,  the  medium  has  degraded  the  apparent  contrast  of  a  high-contrast  (i.e., 
100-percent  modulation)  object  much  below  7  percent  (3.5-percent  modulation) 
at  8  cycles /mrad.  We  do  not  have  data  on  camera  performance  at  very  low 
contrast,  but  one  might  extrapolate  that  at  high  light  levels  one  could  resolve 
a  few  cycles /mrad  at  2 -percent  image  modulation.  This  threshold  level  is 
indicated  in  Fig.  1. 

It  is  important  to  note  that  imaging  systems  have  a  limiting  resolution, 
even  at  high  light  levels.  It  follows  that  increasing  the  illumination  indefi¬ 
nitely  does  not  continue  to  improve  the  resolution. 
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Finally,  if  the  MTF  of  the  medium  is  indeed  as  shown  in  Fig.  1, 
resolution  cannot  be  significantly  improved  by  reduction  of  the  camera  thresh¬ 
old.  This  limitation  is  a  result  of  the  steep  drop  of  the  MTF  shown  in  Fig.  1, 
which  indicates  that  the  image  modulation  below  2  percent  falls  another  10  dB 
at  only  negligibly  higher  spatial  frequencies.  Thus  a  10-dB  improvement  in 
the  limiting  resolution  of  the  camera,  will  not  result  in  significantly  better 
image  quality  through  these  waters. 

The  inherent  resolution  of  the  receiver  is  attained  when  the  wavefront 
from  each  resolvable  element  of  the  object  has  transverse  coherence  and  uni¬ 
form  phase  and  amplitude  across  the  whole  receiver  aperture.  The  coherence 
of  the  radiation  illuminating  the  object  is  important  only  in  the  formation  of 
holograms  or  for  the  special  case  of  observation  of  a  specular  object.  For 
incoherent  imaging  (e.g.,  photography  and  television),  on  the  other  hand, 
coherent  illumination  will  result  in  interference  effects  (e.g.,  speckle  pattern) 
that  can  degrade  the  image.  Such  effects  will  be  observed  if  the  wavefront 
reaching  the  object  exhibits  coherence  over  separations  greater  than  that 
corresponding  to  the  size  of  a  resolvable  element  in  the  receiver.  The  inten¬ 
sity  of  the  illumination  at  the  object  is  a  function  of  the  collimation  or  focusing 
of  the  beam;  in  either  case  the  intensity  is  maximum  when  the  source  exhibits 
transverse  coherence  over  separations  equal  to  its  exit  aperture. 

If  the  medium  dominates,  the  dimensions  of  the  region  of  the  object 

coherently  illuminated  may  be  taken  as  the  e  *  point  of  the  MTF  curves  in 

Fig.  1.  For  the  parameters  plotted,  these  values  are  between  0.05  and  0.5  cm 

(using  the  upper  scale  of  the  abscissa).  In  order  that  interference  phenomena 

not  introduce  additional  noise,  the  angular  resolution  of  the  receiver  should  be 

greater  than  this  value  at  the  object.  At  2  0  m  range  this  corresponds  to  an 

-  5  -  4 

angle  between  3  X  !0  and  3X10  rad. 

B.  BEAM  PATTERNS 

Another  quantity  of  interest  is  the  average  spatial  intensity  distribution 
as  a  function  of  9,  where  the  angle  6  is  indicated  in  Fig.  2. 
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The  average  intensity  at  observation  point  P  of  a  laser  beam 
described  by  an  aperture  function  U(r)  in  the  plane  R  =  0  is  given  by 

<I(P)>  =  e-aR(-^jZ  f  d2£Ms(£,R)  e-(ik/R)£‘£ 

x/U(L+  j£)U*(r  -  -£_)  e(ik/R)£‘£  d2£  (9) 

^The  quantity  a  in  Eq.  (9)  is  the  spectral  volume  attenuation  coefficient 
[the  sum  of  the  volume  absorption  coefficient  and  the  volume  scattering 
(510  rad)  coefficient].  The  effect  of  these  absorptions  and  scatterings  at 
range  R  is  taken  into  account  in  expressions  involving  the  square  of  the  field 
by  the  factor  exp(-a'R).  This  factor  is  to  be  understood  and  is  omitted  in  the 
following.  The  large-angle  scattering  and  absorption  result  in  a  loss  of  power 
from  the  beam.  There  is  no  way  of  distinguishing  between  absorption  and 
large -angle  scattering  (>10  rad)  by  measuring  optical  properties  in  the 
forward  direction .j  As  an  example,  consider  a  plane  wave  of  uniform  ampli¬ 
tude  and  phase  emanating  from  a  circular  aperture  of  diameter  D;  Eq.  (9) 
becomes 


<I(P)>  =  ^  p2  J'  xJQ(2ax)Mg(Dx,  R)Tp(x)  dx 


(10) 


where  Mg  is  the  spherical-wave  MTF  for  the  medium,  x  =  p/D,  P  =  kD  /4R, 
and 


r?(x) 


cos  '  (x) 


for  x  £  1 


d9 


sin  [  2Px(cos  9  -  x)  ]  {l  -  cos[23x(cos  9  -  x|} 


(2  3  x  cos  9 


(23x  cos  9)* 


0  for  x  >  1 


(ID 
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JQ  is  the  zero-order  Bessel  function,  and  cn  =  kDp/2R  -  kD  tan  9/2  =  kD0/2, 
since  0  «  1.  Normalizing  the  intensity  at  a  fixed  range  R  to  the  value  on  the 
axis  (0  =  0),  we  obtain 


<!»(«»  ■ 


g  xJQ(2o'x)Mg(Dx} R)rg (x)  dx 
xrg  (x)Mg(Dx,  R)  dx 


The  Fresnel  and  Frat.ihofer  regions  are  for  values  of  P  greater  than  or  less 
than  one,  respectively.  We  note  that  for  \  -  0.5  pm,  D  -  2  cm  and  R  i  100  m, 
P  £  17.  That  is,  for  all  cases  of  interest  we  are  in  the  Fresnel  region  of  the 
transmitting  aperature. 

The  spherical  wave  MTF  is  given  by  Eq.  (2).  Substituting  this  into 
Eq.  (10)  yields  (see  Appendix) 

<IN(0)>  =  exp(-02/202)  (13) 


where 


02  =  <An  >H 


In  Fig.  3  we  plot  the  normalized  intensity  as  a  function  of  0  for  a  =  50 
cm  and  various  values  of  R  and  <An2>.  The  quantity  0q  is  the  angle  corres¬ 
ponding  to  the  standard  deviation  of  the  gaussian  beam  pattern  and  is  of  the 

order  of  few  tenths  of  a  milliradian,  it  is  plotted  in  Fig.  4  for  various  values 
2 

of  (A  n  )  and  a.  The  curves  presented  here  indicate  the  values  of  beam  spread 

one  w">uld  expect  to  observe  when  imaging  resolution  of  a  few  tenths  of  a 

milliradian  is  achieved.  On  the  other  band,  the  half- power  cone  angle  in  the 

absence  of  index-of- refraction  fluctuations  is  of  the  order  X./D,  i.e.,  a  few 
-  2  -  5 

limes  10  mrad  for  \  ~  5  X  10  cm  and  D  of  some  centimeters.  Hence,  for 

2  -9 

R  ">f  some  tens  of  meters  and  <An  >  £  10'  ,  the  beam  is  spread  by  one  order 
of  magnitude  by  the  large-scale  index-of- refraction  fluctuations. 
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igure  4.  Angle  0Q  as  Function  of  R  for  Various 
Values  of  a 


o  corresponds  to  the  standard  deviation  of  the 
aussian  beam  pattern.) 


From  Fig.  3,  we  note  that  if  propagation  were  determined  by  large-scale 
temperature  and  salinity  fluctuations  of  the  magnitude  assumed  and  if  the 
source  were  perfectly  coherent,  the  half-power  beam  width  even  at  a  range 
of  20  m  would  lie  between  0.06  and  0.2  mrad.  This  corresponds  to  a  coherent 
radiator  in  vacuum  of  between  0.3-  and  0.9-cm  diam.  Therefore,  for  work 
at  20  m,  there  appears  to  be  no  need  for  greater  transverse  coherence,  and 
at  greater  ranges  the  requirement  is  reduced  still  further. 

As  a  possible  experimental  method  of  measuring  the  MTF,  we  note 
that  Eq.  (10)  is  merely  the  Fourier-Bessel  transform  of  the  quantity 
(3  Mg(Dx,  R)Fp(x),  which  can  be  inverted  to  yield 


D2rp(p/D)Ms(p,R)  =  2  it 


<I(p,R)>  p  dp 


(15) 


Hence,  if  the  relative  intensity  <I>  is  known  experimentally  (say,  from  a 
microoensitometer  technique  applied  to  a  direct  exposure  of  a  photographic 
plate  to  the  laser  beam  at  range  R),  one  obtains  directly  the  quantity  Mg 
[by  numerically  integrating  the  righthand  side  of  Eq.  (15)].  The  quantity 
Fp(p/D)  is  essentially  the  MTF  of  the  transmitting  aperture  and  is  assumed 
known  (or  can  be  measured  directly). 

The  MTF  obtained  by  integration  of  Eq.  (15)  can  be  compareu  to  the 
theoretical  expression 


t  <->\  ( P^k2(An^R\ 

Ms(p.R)  =  exp(-* - ^ - j 


if  the  quantities  a  and  ^An  ^  are  measured  simultaneously  with  (1(8)^.  We 
remark  that,  with  6q  defined  by  Eq.  (14),  the  MTF  can  be  written 


Mg(p,  R)  =  exp 


as 


Therefore,  a  measurement  of  the  angle  9  ,  where  the  relative  intensity  is 
down  by  i  A/e  in  comparison  with  its  on-axis  value),  results  in  a  determination 
of  the  MTF.  The  MTF  so  obtained  can  then  be  compared  with  other  simul¬ 
taneous  measurements  of  the  MTF.  Consequently,  these  functions.  Mg  and 
provide  a  means  of  tying  together  independent  measurements  of  resolu¬ 
tion,  beam  divergence,  and  temperature  and  salinity  fluctuations. 

C.  LIMITING  RESOLUTION 

In  conclusion  we  present  an  expression  for  the  limiting  resolution  in 
clear  ocean  water.  The  imaging  properties  of  the  water  are  described  by  the 
MTF.  In  the  model  suggested  here,  the  MTF  is  given  by  Eq.  (2).  The  trans¬ 
verse  distance  p  in  this  expression  is  related  to  the  spatial  frequency  f  by 

f  =  ^  (cycles/rad) 


Consequently,  the  limiting  resolution  is  independent  of  wavelength. 
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APPENDIX.  EVALUATION  OF  AN  INTEGRAL 


We  wish  to  compute 


I-/*  xjQ(2o'x)Mg(Dx,  R)rp(x)  dx 


(A-l) 


where  Fp  is  given  by  Eq.  (li)  and 


Ms(Dx,  R)  =  exp^- 


n2  2.  2/  2vn) 

D  x  k  (An  >R 


(A-2) 


Substituting  Eq.  (A-2)  into  Eq.  (A-l),  we  obtain 


-f,  xJo 


2  2 

(2n'x)rp(x)  exp(-o  x  )  dx 


(A-3) 


where 


Y2  =  k2D2R<An2>/12a 


5  - 1 

For  all  cases  of  interest,  y  »  1  (for  k  ~  10  cm  ,  D  =  2  cm,  R  =  100  m, 
a  ~  50  cm,  and  <An2>  ~  10"9,  Y  ~  105).  Let  y  -  Yx;  then 


I  =  Y 


2/  M^v) 


exp(-y  )Tp(y/Y)  dy 


!«•»  page 
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The  integrand  is  nonzero  only  for  y  less  than  of  the  order  unity.  Since 
V  »  1,  we  may  replace  the  upper  limit  of  the  integral  by  <x>  and  substitute 
Fp(0)  (=  ir/4)  for  Tp(y/Y).  v/e  obtain 


-^2 /  yJ0(T y)  exp<'y2)  dy 


=  exp(-o2/y2) 

8Y 
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